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Metal enrichment of the intracluster medium: SN-driven galactic winds 
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We investigate the role of supernova (SN)-driven galactic winds in the chemical enrichment of the intracluster medium 
(ICM). Such outflows on galactic scales have their origin in huge star forming regions and expel metal enriched material 
out of the galaxies into their surroundings as observed, for example, in the nearby starburst galaxy NGC 253. As massive 
stars in OB-associations explode sequentially, shock waves are driven into the interstellar medium (ISM) of a galaxy and 
merge, forming a superbubble (SB). These SBs expand in a direction perpendicular to the disk plane following the density 
gradient of the ISM. We use the 2D analytical approximation by Kompaneets (1960) to model the expansion of SBs in 
an exponentially stratified ISM. This is modified in order to describe the sequence of SN-explosions as a time-dependent 
process taking into account the main-sequence life-time of the SN-progenitors and using an initial mass function to get 
the number of massive stars per mass interval. The evolution of the bubble in space and time is calculated analytically, 
from which the onset of Rayleigh-Taylor instabilities in the shell can be determined. In its further evolution, the shell will 
break up and high-metallicity gas will be ejected into the halo of the galaxy and even into the ICM. We derive the number 
of stars needed for blow-out depending on the scale height and density of the ambient medium, as well as the fraction of 
alpha- and iron peak elements contained in the hot gas. Finally, the amount of metals injected by Milky Way-type galaxies 
to the ICM is calculated confirming the importance of this enrichment process. 
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1 Introduction 

The space between the galaxies in a cluster is filled with 
a low-density, high-temperature plasma, the so-called intra- 
cluster medium (ICM). It constitutes the largest part of the 
visible mass of a cluster. X-ray observations reveal that the 
ICM contains heavy elements like iron, silicon and oxygen 
instead of existing entirely of primordial gas. On average, 
metallicities range from 1/3 to 1/2 of the solar abundances 
(e.g. Renzini 1997, Fukazawa et al. 1998) (see Fig.Q). Al- 
ready in 1975, Larson and Dinerstein discussed the chemi- 
cal enrichment of the ICM in their theoretical paper and sug- 
gested that mass loss from galaxies in a cluster increases the 
metallicity of the ICM. Shortly after that, the iron emission 
line was detected in the X-ray spectra of the Perseus cluster 
(Mitchell 1975). A large number of samples of nearby clus- 
ters (e.g. Mushotzky et al. 1978; White 2000) has been ana- 
lyzed until now, as well as clusters out to redshifts of z ~ 1 
(e.g. Tozzi et al. 2003, Balestra et al. 2007). Only a small 
evolutionary trend can be seen in the abundances. It is obvi- 
ous that investigations of the enrichment of the ICM help us 
to understand the origin of the ICM and the star formation 
history of the universe. Since heavy elements are the prod- 
uct of nucleosynthesis inside the cluster galaxies, the ques- 
tion comes up, how they can be transported from the cluster 
galaxies into the ICM. Possible enrichment processes are 
galactic winds, ram pressure stripping and jets from active 
galactic nuclei. Exploring galaxy clusters with X-ray obser- 
vatories like XMM-Newton results in 2D abundance maps, 



which show, in general, an inhomogeneous distribution of 
the metals (Durret et al. 2005). Studying the distribution 
of metals in the ICM provides constraints on the efficiency 
of these processes and on the evolution of clusters and the 
galaxies within (Schindler et al. 2005). We need to know 
not only how these processes are working but also must cal- 
culate the amount of metals that are ejected or stripped in 
order to learn about the formation and chemical evolution 
of galaxies. 

2 Galactic outflows 

Galactic winds are powerful outflows of galaxies originat- 
ing from massive bursts of star formation. Thus they are a 
common phenomenon in starburst galaxies, where the star 
formation rate is enhanced by a factor of 10 or more. In 
addition cosmic rays can assist in driving galactic winds, 
provided that the coupling between high energy particles 
and gas is sufficiently strong (Breitschwerdt et al. 1991, Ev- 
erett et al. 2008). Heavy elements are synthesized in mas- 
sive stars, which end their life as supernovae. The released 
energy of successive SNe in star forming regions drives the 
outflow and the metal enriched material is expelled into the 
halo of a galaxy or even into intracluster space. 
Galactic winds as an enrichment process were proposed by 
Larson & Dinerstein (1975) and then by De Young (1978). 
The amount of mass loss depends on properties of a galaxy 
itself, such as mass and disk scale length, and on the en- 
vironment. An important question is, where the transition 
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Fig. 1 Iron abundances as a function of ICM temperature 
for clusters and groups. Clusters at moderately high redshift 
(z) ~ 0.35 are shown by small filled circles (Renzini 1997). 



of metal ejection to metal retention is. Wiebe et al. (1999) 
find that galaxies with more than 5 • 10 11 Mq do not eject 
any material at all. Dwarf galaxies may suffer a complete 
blow-away due to their shallow potentials, while galaxies 
with 10 9 M Q lose less than 3 % of the metal enriched gas 
(MacLow & Ferrara 1999). With outflow speeds of about 
200-1500 km/s (Bland-Hawthorn, Veilleux k Cecil 2007), 
the escape velocity from the gravitational potential is easily 
exceeded in low-mass galaxies, such that the gas is lost from 
the galaxies. The other possible fate of the high-metallicity 
gas is that it cools and falls back onto the galaxy, becom- 
ing part of the galactic fountain (Shapiro & Field 1976) and 
consequently returning the newly synthesized material to 
the ISM. Contrary to ram pressure stripping, galactic winds 
are usually suppressed in the central regions of a cluster due 
to the high ICM pressure, against which outflows have to 
work. 



3 Superbubble dynamics 

Highly energetic interaction processes between the compo- 
nents of the ISM are found where strong explosions like 
supernovae occur or stellar winds from massive stars sweep 
up the ISM. Massive stars are usually found in associations 
rather than randomly in the field (Lada & Lada 2003). They 
are strongly correlated in space in regions of ~ 100 pc radius 
and in time. The stellar winds from these stars and multi- 
ple SN-explosions produce shock waves and carve out large 
holes filled with hot, tenuous and metal enriched gas. These 
SBs can blow out of the disk and expel the hot gas into 
the halo. Moreover, they deliver energy to the ISM and into 
the halo, and hence provide important information about the 
heating of the disk and the galactic corona by SNe. 
The aim of our work was to develop a time-dependent mo- 
del for superbubble evolution to find out if galactic winds 
can be a dominant enrichment process, and to estimate how 
many metals are released in a wind. 



3.1 Kompaneets' approximation 

At early times in the evolution of a superbubble the sur- 
rounding medium is assumed to appear homogeneous and 
therefore the surface of the shock front remains spherical. 
The similarity solution (Sedov 1946, Taylor 1950) can be 
used to describe the expansion of the bubble in these early 
stages as long as the radius is less than the ambient den- 
sity scale height. Due to the negative density gradient the 
bubble becomes elongated after reaching one scale height 
above the galactic plane (Fig. 0. We use the 2D analyti- 
cal approximation found by Kompaneets (1960), which - 
in its original form - describes the propagation of an adia- 
batic shock wave generated by a point explosion in an ex- 
ponentially stratified medium. The density variation is given 
by p{z) = po exp(— z/H). Here, z is the height above the 
galactic plane, po is the density at z — and H is the scale 
height of the interstellar gas. In our work, we use Kompa- 
neets' approximation (hereafter KA) to determine the shape 
and the time evolution of a superbubble analytically, which 
provides a proper insight into the physics of bubble evolu- 
tion. 

The axially symmetric problem is described in cylindrical 
coordinates (r, z). In order to use the KA for the investiga- 
tion of SB evolution the following assumptions have to be 
made (see also Schiano 1985; Bisnovatyi-Kogan & Silich 
1995; Basu et al. 1999): (i) a strong shock is propagating 
into the undisturbed ISM. This means that the pressure in 
front of the shock, i.e. of the ambient gas is negligible; (ii) 
the pressure of the shocked gas is spatially uniform and (iii) 
almost all the swept-up gas behind the shock front is lo- 
cated in a thin shell. In the analysis of the evolution of the 
shock surface, the Rankine-Hugoniot jump conditions have 
to be solved at every point of the shock surface (Bisnovatyi- 
Kogan & Silich 1995; Basu et al. 1999). The solution is 
performed using a transformed time variable (in units of a 
length), which represents the evolutionary status of a bubble 
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with E t h being the thermal energy of the shocked gas and 
£1 being the volume confined by the shock front. A ratio of 
specific heats 7 = 5/3 is used. The result is the half-width 
extension of the bubble parallel to the galactic plane (see 
Fig.0 
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Top and bottom of the bubble are given by 

z u {y) = -2H ■ In (1 - y/2) 

and 

z d (y) = -2H-ln(l+y/2), 
respectively, with y = y/H. 
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Fig. 2 Position of the shock front at y = 0.5, 1.0, 1.4, 1.7, 
1.9, 1.98, 1.99, and 2.0. 



3.2 Superbubbles in stratified disks 

In order to derive analytical expressions for time, velocity 
and acceleration of the superbubble we approximate the vol- 
ume of the bubble by an ellipsoicQ, with a semimajor axis 
a{y) = (z u {y) - z d {y))/2 and a semiminor axis r max (y) = 
2H arcsin (y/2). Thus, we get the bubble volume V(y) = 
4^ a(y) r^ ax as a function of the time variable. The devia- 
tion from the numerical integration of the volume amounts 
to only ~ 2 % at very late stages of evolution (y = 1.9). 
In our description of superbubble evolution we want the ini- 
tial point-like energy deposition in the KA to be replaced 
by a time-dependent energy input rate. We calculate the en- 
ergy thermalized at the inner shock following Weaver et 
al. (1977), but using the time-dependent energy input rate 
LsB(t) = Limf ■ t S instead of a constant wind. In such a 
model, masses and numbers of stars are given by an Initial 
Mass Function (e.g. Fuchs et al. 2006) 



dN{m) 
dm 



N ■ m?. 



(5) 



With r = 7 + 1 and by fixing a stellar mass range, the 
integration of the equation above results in the total number 



1 Also Maciejewski and Cox (1999) approximate the shape of the shock 
front by a prolate ellipsoid. 
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Fig. 3 Dimensionless velocity (solid line) at the top of the 
bubble function of y. The minimum of the velocity 

is reached where y acc = 0.68 (dashed line). 

of stars between a lower mass limit mi and an upper mass 
limit m u 

No n 

777, 1 



with masses in solar units. The normalization factor Nq is 
given by the number of stars belonging to the OB-associa- 
tion. We obtain it by fixing the number of stars in the last 
mass bin (m„ — 1, m u ) to be Aob = 1- 
Treating the number of stars between (m, m + dm) as a 
function of mass and expressing the stellar mass by its main 
sequence lifetime (i, t—dt) in order to obtain a time-sequen- 
ce of exploding stars (Berghofer & Breitschwerdt 2002) re- 
sults in 



L 3 B(t) = E. 
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The main sequence life time for massive stars as a function 
of mass (in units of seconds) can be described by a power- 
law t(m) = Km~ n . Values taken from Fuchs et al. (2006) 
are k = 1.6 • 10 8 yr and rj = 0.932. Rearranging leads to 
the mass of a star expressed by its main sequence lifetime 
m(t) = {t/n)- 1 ^. We get 

isB(i) = ^^(£y' ^+,, . (8) 

77 • K \ K J 

Since Lsb(£) ^ s k we nave ^ = ~ (§ + !)• Finally, the 
result for the thermal energy as a function of time is 
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75 + 11 

Checking for a constant energy input rate with 5 = corre- 
sponds to the Weaver et al. (1977) solution with a constant 
factor of 5/11. We use the energy input rate, the density and 
the scale height of the interstellar medium to define a time 
scale in units of seconds t = (p ■ i/ 5 /Li M F) (1/(<5+3)) . By 
rearranging equ.Q]we obtain the time as a function of y 

'{s + 3) 2 {7d + ny " ' 
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Fig. 4 Instability time scale (solid line) and dynamical 
time scale (dotted line) for n = 0. 1 cm -3 , H = 500 pc, Nob 
= 20 as a function of the time variable y. Acceleration starts 
at y acc — 0.68, the curves cross each other at y r t; = 0.70 
(indicated by the dashed line). 
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Fig. 5 Fragmentation time scale at the top of the bubble 
as a function of the total number of stars for associations 
with 10-3000 OB-stars, again for ISM-values of the Lock- 
man layer. Dashed lines, designating certain time intervals, 
are used for illustrative purposes. 



with a = 1/(8 + 3) and /3 = J ■ By expressing N as 
a function of Nqb an d the integral of the square root of the 
volume by a simple power law, we derive the time in typical 
units of the ISM 
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The errors for this approximation are below 2% (y — 1.9, 
~ 100 stars, Lockman layer). 

Velocity z u (y) and acceleration z u (y) of the top of the bub- 
ble are derived by calculating the first and second time deri- 
vative, respectively, of the coordinate z u (equ.0. Using dif- 
ferent formulas for the main sequence life time or also vary- 
ing the slope of the IMF may change the value of y acc and 
hence the point in time where the acceleration sets in. For 
an IMF with slope Y = -1.35 (Salpeter 1955) we see that 
the velocity has its minimum at j/ acc = 0.68 (Fig. [3]). By 
assuming that the shock has to stay strong all the time we 
introduce our blow-out condition: the velocity of the outer 
shock needs to have a Mach number of at least M > 5. We 
find that blow-out into the Lockman-layer (Lockman 1984) 
of a Milky Way-like galaxy (n = 0.1 cm -3 , H = 500 pc, and 
T = 6000 K) requires an association with about 20 OB-stars, 
which, incidentally corresponds to the initial stellar content 
of the Local Bubble. 

The acceleration leads to the development of Rayleigh-Tay- 
lor instabilities (RTIs) at the contact discontinuity, since a 
dense shell accelerated by the tenuous gas inside the bubble 
constitutes an unstable hydrodynamic configuration. RTIs 
grow exponentially on a characteristic timescale t It i, which 
is calculated for the top of the bubble (in units of seconds) 



*rti,z u (y) 



I Ad e (y) _ 4p e-^(v) + Pin (y) 
2Trz u (y) ' 4p e-^(y) - p in (y)' 



(12) 



The wavelength of the perturbation is comparable to the 
shell thickness Ad e and p ln means the density of the hot 
shocked gas. 

If the timescale of an instability to grow is smaller than the 
dynamical time-scale idyn(y) = o,(y)/z u (y) at any point in 
time, the shell starts to break up at this moment. When com- 
paring these two timescales, the value of y can be found, 
where t r ti, Zu < ^dyn- The break up of the shell starts for 
all bubble configurations shortly after acceleration sets in, 
which happens at y Itl — 0.7 for rich and poor associations. 
Fig.|4]shows the behavior of these two curves for an associ- 
ation with 20 OB-stars. In our calculations it is assumed that 
the exponentially growing instability is fully developed af- 
ter approximately 3 • t T tu which we call fragmentation time 
scale (Fig. 0. At this point, the instability has reached the 
fully non-linear regime. The shell will break-up and the hot 
material inside the bubble can escape and enrich the sur- 
rounding medium with heavy elements. We have to take 
into account that fragmentation has to occur in a reasonable 
time scale, before or not long after all SNe have exploded 
and where shear motions or turbulence in the ISM have no 
great influence on the shell structure. If we define such a 
time interval to be 30 Myr after the acceleration at top of 
the bubble had started, an association needs to have at least 
32 OB-stars. In that case, the shell is fragmented at z u about 
41 Myr after the first SN-explosion happened. 



4 ICM metallicities 

In this section, we calculate the amount of metals ejected 
in a wind, in order to compare it to the metal content of a 
galaxy cluster. Yields for the elements iron ( 56 Fc) and oxy- 
gen ( 16 0) from Woosley & Weaver (1995) are adopted to 
a stellar mass range between 8 and 120 M®. For the iron 
mass produced by a massive star with mass m we find a 
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Fig. 6 Iron mass ejected by SNe in OB-associations. The 
mass is shown as a function of the total number of OB-stars 
in the star cluster using an IMF with V = —1.35. 

power law (given in solar masses) 

M Fe (m) = 0.02 • m ' 27 - 0.03. (13) 

A similar result is obtained for oxygen 

Mo(m) = 0.049 • m 1162 . (14) 

Integration over the whole mass interval of an OB-associa- 
tion yields the total mass of iron and oxygen, respectively, 
contained in the hot interior of a SB. Fig. |6]shows the iron 
mass ejected by SNe as a function of the total number of 
OB-stars, Fig.|7]shows the relation for oxygen. 
Now, that we know the mass of iron and oxygen for SBs 
with different richness, the next step is to consider the total 
fraction of metals expelled by SNe in a whole galaxy over 
its evolution. 

Assuming a SFR of 5 M©/yr gives - over an interval of 
10 10 yr - about 5 • 10 10 M of material which is used to 
form stars with masses between 0.08 and 120 Mq. Addi- 
tionally, about 10 11 M are supplied by 5 strong starburst 
episodes (2OOAf /yr) where each episode lasts 10 s yr. Al- 
together, there are 1.5 ■ 10 11 M representing a Milky-Way 
like galaxy. The fraction of mass found in SN-progenitors is 
about 13 %. The available mass of Mob = 1.9- 1O 1O M for 
the formation of high-mass stars has to be distributed among 
galactic OB-associations. Observations show that the num- 
ber of OB-associations in spiral galaxies follows an Initial 
Cluster Mass Function (ICMF) with a slope U, = -0.75 
(Dowell, Buckalew & Tan 2008). The number of associ- 
ations as a function of their total mass can be calculated 
now. High mass star formation is assumed to happen only in 
these associations. The total mass Mob of massive stars is 
distributed among associations that do not have more than 
3000 OB-stars or a mass M u = 6.4 ■ 10 4 M in form of 
OB-stars. The least massive association has two OB-stars 
or Mi ~ 20 M . 

First, a normalization constant has to be calculated 
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Fig. 7 Same as Fig. |6l but for oxygen. 



Then the total number of OB-associations can be derived in 
this mass interval by integration over the ICMF 

AWob(M) = • (M 2 n - Mf) . (16) 

For Mi = Mi and M 2 = M u , we get a total number of 
OB-associations of ~ 4.8 • 10 7 . Since only SBs with at least 
32 OB-stars or M\ = 400 M are taken into account (de- 
rived from the fragmentation timescale), the number of OB- 
associations that may enrich the surrounding intergalactic 
medium is ~ 5 ■ 10 6 . If we sum up the number of OB- 
associations with a certain mass multiplied with the iron or 
oxygen mass per cluster, we get the amount of heavy ele- 
ments that is lost to the galactic neighborhood. We find that 
1.46 • 10 7 Mq of iron and 1.3 • 10 9 M of oxygen produced 
by massive stars are expelled by a Milky Way-type galaxy 
over 10 10 yr. 

However, a significant contribution to the blown-out metals 
comes from the entrainment of ISM, which has been pol- 
luted by generations of previous star clusters. Thus we need 
to calculate the metal amount in the swept-up shell of the SB 
as well. We estimate the abundances in the ISM and thus in 
the shell to have, on average, about half of the solar metallic- 
ity over the entire evolution of the Galaxy (e.g. Ibukiyama 
2004). We get the mass of the shell of the swept-up ISM 
by integrating the ellipsoidal shape over the density gradi- 
ent from Zd to z u . At the moment of onset of RTIs, there 
are typically ~ 5 ■ 10 5 M of gas in the shell. We derive 
from geometrical arguments that the mass of the shell frag- 
ments that are lost in the outflow, i.e. the cap of the bubble, 
amounts to about 10 % of the total shell mass. Using abun- 
dances by Grevesse, Asplund & Sauval (2007), we derive 
that 30 Mq of iron and 135 Mq of oxygen per SB are in- 
jected to the ICM that way. Altogether, the shell fragments 
of one galaxy yield an iron mass of about 1.5 • 10 8 M and 
an oxygen mass of about 6.75 • 10 8 Mq. 
Our intention is to compare these results to a cluster of 100 
Milky Way-like galaxies. Such a cluster has a total visible 
mass of about 1.2 ■ 1O 14 M . With 1/3 solar metallicity, 
there are ~ 4.4 . 10 10 M 0i f c in a galaxy cluster. Oxygen 
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abundances are 2/3 of the solar abundances, thus the oxy- 
gen mass in the cluster is found to be ~ 4.1 • 10 11 M 0) o. 
Metals in the shell and in the ejecta yield, finally, for 100 
galaxies 1.65 • 10 10 M 0iFe and 1.98 • 10 11 M 0! o- Compar- 
ison of our results to the estimates derived from observa- 
tions, shows that more than 1/3 of the iron and almost 50 % 
of the oxygen mass in the ICM can be explained by galactic 
outflows. In both cases, metals in the fragmenting shell as 
well as those ejected by SNe contribute to the enrichment of 
the intracluster gas. 

5 Closing remarks 

1 . Although the KA does not take into account inertia of 
the thin shell, cooling of the shocked gas and gravity, a 
good agreement between the bubble geometry derived 
from this model and that of observed SBs usually is 
found (e.g. Gulliford 1974, Basu et al. 1999). The ap- 
proximation seems to be a good choice to describe SB 
evolution in a simple, but realistic way, but instead of a 
pure exponentially decreasing density law one symmet- 
ric to the midplane should be included. 

2. Comparison to other IMFs (e.g. T = -1.15 and T = 
— 1.7) reveals that a similar number of stars is needed 
for blow-out into the Lockman layer (18 and 26 stars, 
respectively). 

3. The calculations of the ICM metallicities, although ba- 
sed on a rough estimate only, show convincingly that 
SN-driven galactic winds are a relevant source of metal- 
enriched material and play an important role in the en- 
richment of the ICM. 
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